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Micro heat engines have attracted considerable interest in recent years for their potential exploitation as 
micro power sources in microsystems and portable devices. Thermodynamic modeling can predict the 
theoretical performance that can be potentially achieved by micro heat engine designs. An appropriate 
model can not only provide key information at the design stage but also indicate the potential room for 
improvement in existing micro heat engines. However, there are few models reported to date which are 
suitable for evaluating the power performance of micro heat engines. This paper presents a new ther¬ 
modynamic model for determining the theoretical limit of power performance of micro heat engines 
with consideration to finite heat input and heat leakage. By matching the model components to those of 
a representative heat engine layout, the theoretical power, power density, and thermal efficiency 
achievable for a micro heat engine can be obtained for a given set of design parameters. The effects of key 
design parameters such as length and thermal conductivity of the engine material on these theoretical 
outputs are also investigated. Possible trade-offs among these performance objectives are discussed. 
Performance results derived from the developed model are compared with those of a working micro heat 
engine (P3) as an example. 

© 2011 Elsevier Ltd. All rights reserved. 


1. Introduction 

There is an ever increasing demand for power sources on the 
micron scale which are essential for powering a wide range of 
microsensors, microactuators, and other applications [1], In 
contrast to macro scale power sources, micro power supplies are 
compact, lightweight, and promise to deliver high power density 
[2], Thus, in recent years, there has been significant research 
worldwide in the development of a variety of micro power gener¬ 
ation systems [3,4]. Among these systems, micro heat engines 
attract major interest as they can be fabricated using mainstream 
MEMS processes and thus potentially integrated with a wide 
variety of microsystems [4], In addition, the combustion of fuel in 
micro heat engines results in high energy density power sources 
even when heat engine conversion losses are taken into account 
[5], Typical micro heat engine development reported to date 
includes the MIT micro gas turbine [6], the P3 micro engine [7], 
micro Wankel engine [8], micro steam engine [9], and micro 
reciprocating engines [10]. 
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Similar to that of their macro scale counterparts, the overall 
performance of micro heat engines is determined by tradeoffs of 
a number of parameters which need to be considered at the design 
stage. Thermodynamic modeling has been used to predict theo¬ 
retical limits, e.g. Carnot limit of maximum thermal efficiency in 
general heat engines. While it is unlikely that real micro heat 
engines will operate at the theoretical limits predicted by ther¬ 
modynamics, evaluation of the thermodynamic limit is important 
as a measure of the potential room for improvement in real-life 
micro heat engines. Furthermore, theoretical models can also be 
used to evaluate potential performance with different design 
parameters. In particular, the down-scaling effect for regenerative 
heat engines in terms of thermal efficiency has been reported [11]. 
However, to the best of our knowledge, there are no systematic 
modeling studies on the theoretical limit of power and power 
density of micro heat engines as functions of length scale and 
material selection. Moreover, as thermal isolation of the hot and 
cold sections of an engine becomes increasingly challenging at the 
micron scale [12], heat leakage is thus one of the key considerations 
in micro engine designs [13]. Previous studies in miniaturized heat 
engines have been based on a thermodynamic heat engine model 
with leakage [14], as shown in Fig. 1(a), where the relevant 
parameters are listed in Table 1. Although the thermal efficiency vs. 
leakage characteristic of a micro heat engine was derived from this 
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Fig. 1. (a) Thermodynamic model of heat engine with leakage (b) New finite heat 
model of micro heat engine. 

model, the effect of leakage on power generation was not predicted 
[11]. This is because an infinite amount of input heat Qh was 
assumed to be drawn from the source Th, leading to generated 
power being independent from leakage [15]. However, in real heat 
engines, there is a limit in the heat flux that can be drawn from 
a source [16], which is particularly important for micro heat 
engines due to the down-scaling of the cross-sectional area and 
thus much reduced heat flow in total. This can play a critical role in 
determining the power and power density output of micro heat 
engines which was not considered in previous models. 

This paper focuses on the development of a new thermody¬ 
namic model which includes both the leakage conductance and an 
additional source conductance (Ks) between the source and the 
engine emphasizing the finite heat input in real situations. Devel¬ 
opment of this model is illustrated by matching the model 
components to those of a representative heat engine layout. This 
allows the evaluation of the theoretical limits of power, power 
density, and thermal efficiency of micro heat engines as well as 
their dependence on key design parameters which can be poten¬ 
tially exploited for performance optimization of micro heat 
engines. 

Table 1 

Nomenclature. 

Th 
Tc 
Ts 
Tin 
T 1C 
Ks 
Kh 
K c 
Kl 
Qs 
Qh 
Qc 
Qi 
Or 
W 

T 

k L 

l<s 

L 
A 
to 

Ao 
s 


hot side temperature 

cold side temperature 

source temperature 

internal hot side temperature 

internal cold side temperature 

source thermal conductance 

engine hot side thermal conductance 

engine cold side thermal conductance 

leakage thermal conductance 

rate of heat input to system (engine and leakage) 

rate of heat input to engine 

rate of heat rejection from engine 

rate of heat leakage 

total rate of heat input 

rate of work (power) 

ratio of internal temperatures 

material thermal conductivity of leakage path 

material thermal conductivity of source 

length of micro heat engine 

cross-sectional area of micro heat engine 

additional length due to the source conductance 

additional area due to the leakage conductance 

length scale 


2. Finite heat model 

Fig. 1(b) shows a proposed thermodynamic model for micro 
heat engines. It is composed of a source conductance I<s, a leakage 
component I< L , and an endoreversible engine component E, as 
shown inside the dashed box. The engine component E is in turn 
composed of a reversible Carnot engine running between a hot side 
temperature of T !H and a cold side temperature of T,c, with finite 
thermal conductance Kh connecting T,h to the external hot side 
temperature Th, and thermal conductance I<c similarly connecting 
Tjc to the external cold side temperature To The leakage component 
is a thermal conductance I< L running between Tj and To shunting 
heat away from the engine component. The source conductance Ks 
connects the source to the rest of the system, representing a finite 
heat input rate. The variables Qs • Qh • Qc • Ql are cycle-time- 
averaged values of the heat flow rates through the conductances Ks, 
Kh, Kc, and I< L respectively. 

The new component I<s for micro heat engines can be justified as 
follows. Since heat leakage is present in all engines whether large 
or small, it is necessary to include a heat leakage component I< L for 
a realistic engine model. With micro engines however, a further 
issue is that heat input into the engine becomes constricted as 
engine size is reduced, while leakage tends to increase. This causes 
Ki and Ks to become comparable in magnitude in micro heat 
engines. Whereas in a conventional engine, the magnitude of Ks 
would be much greater than that of I<l, and hence Ks can be ignored 
in the thermodynamic model. This argument is further demon¬ 
strated by considering a simplified engine layout as well as a few 
real micro engines. 


2.3. Illustration of I<l/I<s ratio 


Heat losses in real micro heat engines can be from a range of 
components [ 14,17], especially for engines with complex geometries. 
However, in order to illustrate the increasing I< L /K S ratio, we consider 
a simplified layout of an engine chamber as shown in Fig. 2(a), to 
represent Ks and I<l. Note that this layout directly corresponds to the 
model shown in Fig. 1(b). Heat input comes externally from the top 
through the plate I<s and is conducted to both the engine E and the 
sidewalls K L , representing the parallel heat loss through the walls of 
the engine. While the component E is an endoreversible thermody¬ 
namic model that has been used to model both internal and external 
combustion engines, including gas turbines [18], Stirling engines 
[19], and other engine types [16], the I<s component most accurately 
models heat transfer from an external source (i.e. an external 
combustion type engine). As will be discussed later, I<s for other 
engine types may be modeled differently. 

The magnitude of the heat conductances mainly depends on the 
thermal conductivity of the engine construction material, the 
engine cross-sectional area and the length. Using a linear model of 
heat conduction, the conductances I<s and Kl can be calculated as 
follows: 


3Cc = k, 


A + Aq 
; Lo 


(1) 


K l = k L ^ ( 2 ) 

Where the parameters ks and I<l are the thermal conductivities of the 
source and leakage sections of the heat engine, respectively. A 0 is the 
cross-sectional area of the leakage path, while Lo is the thickness of 
the source conductance path. A and I are the overall area and length 
of the heat engine. One way to investigate the effect of scaling is 
through an abstract length scale parameter s which represents size 
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Fig. 2. (a) Illustration of thermal conductances in a micro heat engine (b) Effect of length scale on l( L /K s ratio. 


in one dimension. Thus, length L is proportional to s and area A is 
proportional to s 2 . From Eqs. (1) and (2), it can be seen that down- 
scaling of micro heat engines lead to much reduced ( K s xA<xs 2 ) 
and increased Ki ( I< l xL^ ). In reality, Ao would also somewhat 
scale with s (e.g. if the walls of the engine were kept to constant 
thickness), so for the example layout this leads to a ratio I<lII<s which 
is proportional to s~ 2 as shown in Fig. 2(b). The increasing KJI<s ratio 
with the down-scaling of micro heat engines suggests that heat 
leakage and finite heat input effect becomes relatively more prom¬ 
inent at the micro scale in comparison with that of macro engines 
and thus should be taken into consideration for modeling, which is 
underlined by the developed new model in this work. 

2.2. Applicability to real micro engines 

The P3 engine is a frequently reported [7] [20] micro heat engine 
that features a thermal switch to control heat addition and rejec¬ 
tion. The P3 engine is basically a chamber filled with a phase change 
fluid, one wall of which is a flexible membrane that acts as the 
generator. Heat input into the engine is controlled by a thermal 
switch. Heat loss occurs through the walls of the engine as well as 
through the working fluid. 

In the case of the P3, I<s is clearly the thermal conductance of the 
thermal switch, i.e. the inverse of the thermal resistance (3.8 K/W) 
of the switch in its “on” state [20], i.e. I<$ = 0.263 W/K. Clearly I<$ for 
the P3 has the same area dependence as in Eq. (1 ). 1<l on the other 
hand can be calculated from the geometry and material properties 
using Eq. (2). The reported [2] P3 engine has cross-sectional area 
A = 19.6 mm 2 , enclosing a cavity that is 75 pm thick with sidewalls 
that are composed of layers of semiconductor tape. A working fluid, 
3M PF-5060DL with thermal conductivity 0.057 W/m-K [21] is 
enclosed within the cavity. Due to the relatively small heat transfer 
area of the thermal switch (4 mm square), the heat leakage is 
dominated by heat transfer through the working fluid. The value of 
Ki can thus be calculated from Eq. (2) using the cross-sectional area 
and conductivity of the working fluid, the value of which is found to 
be 1< L = 0.015 W/K, which leads to I< L /I< S = 0.057. 

In comparison, the I< L /I<s ratio of a conventionally sized engine 
such as the GPU-3 can be roughly estimated from its published 
specifications [22]. Ki can be approximated using Eq. (2), from the 
regenerator area of about 0.01 mm 2 per cylinder and length of 
22.6 mm. Kg is in turn approximated using Eq. (1 ), from the heater 
area of about 286.5 mm 2 per cylinder and length of 245.3 mm. Both 


heater and regenerator are made of stainless steel, with thermal 
conductivity 15 W/m-K. This leads to a I( L /K S ratio of 0.00038, 
which is much smaller than that of the P3 example. Hence, finite 
heat input is not a strong consideration for the GPU-3 engine and 
the new model’s prediction of its performance would not differ 
greatly from previous models. 

While not quantified, I<h and l<c in Fig. 3 represent the effective 
heat transfer coefficient into and out of the working fluid, respec¬ 
tively. These values dictate the absolute level of power that the 
engine is capable of generating, as in the conventional model of 
Fig. 1(a). The relationship of I<s and Kl on the other hand, which is 
the focus of this study, represents a penalty to the generated power 
of an engine as it is miniaturized. 

For additional illustration, we consider another engine type, the 
MIT micro turbine [4[. K$ can be interpreted as a combustion 
effectiveness factor. If this micro turbine were of external 
combustion type, then I<$ would simply be the conductance of the 
input heat exchanger and can represented by Eq. (1). Since it is in 
fact powered by internal combustion, the geometric dependence of 



Fig. 3. Power vs engine length at different thermal conductivities, Th = 400K, 
T c = 300K, P max calculated as power generated without leakage 
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Ks for this engine would be more complicated. For example, is 
influenced by parameters such as the chamber size, quenching 
distance, and fluid flow of fuel, all of which tends to reduce effec¬ 
tiveness of combustion as engine size is reduced [1,5]. Ki again is 
the heat loss through the engine walls, and can be approximated 
using Eq. (2), the net effect of which is that the K L /I< S ratio tends to 
increase in the micro scale, which is the same result demonstrated 
in previous illustrations and hence justifies the use of the new 
thermodynamic model for micro engines. 

Clearly, the developed model is most suitable for representing 
external combustion type engines, but can also apply to internal 
combustion engine types with suitable modification of Eqs. (1) and (2). 


Eq. (12) is an expression of power output W as a function of the 
temperature ratio t. The maximum of a function occurs at the point 
when its derivative is zero. Hence, by applying the condition 
9VV 

— = 0, we can obtain the maximum power that can be generated 

01 

by a micro heat engine with certain leakage and finite heat input: 


W m ax 


K h K c 


( VKsTh+KJc - s/KsTc + KJc) 2 
I< s K h + K s Kc + K h K l + K C K L + K h K c 


(13) 


with the following expression for r corresponding to the maximum 
power point: 


2.3. Model equations 

Power and power density are key performance considerations 
for micro heat engines. Based on the new thermodynamic model, 
the theoretical limit of these performance objectives can be eval¬ 
uated with consideration to leakage as follows: 

The heat flow rate Q H represents the part of heat input which is 
usefully converted to work, while Q L represents the leakage that 
dissipates across Ki. Qs is the total input heat rate from the source, 
and is the sum of Q H and Q L . 


Qs = Qh + Ql 

From Fourier’s law of heat conduction: 

(3) 

Q s = K s (T h - T s ) 

(4) 

Q h = K h (T s - T iH ) 

(5) 


Ql = K L (T S - T c ) (6) 

Substituting (4), (5), and (6) in (3), we obtain: 
t ( K s + Kh+Kl)Ts-K s Th-K l Tc 

T[h =- Kh - (7) 

Since the internal Carnot engine is reversible, the ratio t between 
the internal operating temperatures, Tic and T;h can be calculated 
as: 


T 


(K s + K l )T c 
K s T h +K l Tc 


The relevant power density can thus be obtained as: 


(14) 


W W 

Power density = -,—,-= — A --— - —-— (15) 

total volume (A+A 0 )(L + L 0 ) 

By using Eqs. (13) and (15), the effect of the leakage length and 
thermal conductivity on the limit of maximum power and power 
density of micro heat engines can be obtained, which can then 
provide practical design guidance. For instance, since the 
maximum power is dependent on various thermal conductances 
within the heat engine, it suggests that the development of micro 
scale heat exchangers such as the P3 thermal switch [20] could 
have a large effect on the power output. It is important to note that 
in extreme and ideal cases when K s = °° or when I< L = 0 , Eq. (14) 
is reduced to i = y/T c /T H , which is the maximum power condition 
predicted by using the previously reported thermodynamic model 
in which infinite heat input was implicitly assumed [23]. 

Another quantity of interest is the thermal efficiency at 
maximum power q, which can be expressed as the maximum 
power W max over the total heat input rate Q s : 


W _ W 
Pis ~ k s(T h -T s ) 


(16) 


With the value of W max obtained from (13) and Ts from (10). 


T _ T ic _Qc _ jCcOic ~ T c ) 

Tm Qh K h (T s ~ T^) 

Which leads to the expression for T !H as: 


_ tK H T s + K C T C 
iH ~ t(I< h + K c ) 

Joining Eqs. (7) and (9), Ts can be obtained as: 


(9) 


t(K h + K c )(K s T h + K l T c ) + K H K C T C 
t(K s K h + K S K C + K h K l + K c I< l + K h I<c ) 


( 10 ) 


The power output can be obtained from thermal efficiency and heat 
input as: 


w = (1 - T )Q h = (1 - t )K h (T s - T m ) (11) 

by inputting Ts from (10) and Tm from (9), we obtain: 


W = K h K c 


K s T h + 2K l T c + I< S T C - t (I< s T h + K L T C ) - V s + I< L )T C 
K s K h + K S K C + K h I< l + I< c K l + K h K c 


( 12 ) 


3. Analysis and discussion 

3.1. Power 

Since power generation is one of the key applications for micro 
heat engines, it is important to predict and optimize the upper limit 
of maximum power that can be potentially produced. Design 
parameters such as thermal conductivity of the engine construction 
material, length scale, and hot side temperature can all play 
important roles, as indicated by Eq. (13). Fig. 3 shows the power 
output as a function of engine length L. It is observed that for larger 
length scales, the power output is close to P max which represents 
power generated without any leakage. However, with the decrease 
of L, the power output rapidly falls off due to the increase of leakage 
and I< l /K s ratio. The minimum L for a viable device depends largely 
on the thermal conductivity k L . For example, to avoid excessive 
leakage, silicon (ki = 150 W/m-K) is only suitable for engines with L 
in the range of 10—1000 mm; whereas for smaller engines with L in 
the range of 0.1-1 mm, polymers such as PMMA and SU-8 
(ki = 0.1—0.2 W/m-K) are better choices as construction materials. 
Fig. 4 (a) shows normalized power vs. length at different source 
temperatures (Th) where output power increases at higher source 
temperatures. If the micro heat engine is designed for an output 
















2690 


K. Khu et at / Energy 36 (2011) 2686-2692 



Fig. 4. (a) Power vs engine length at different source temperatures, P max calculated as power generated without leakage, (b) L min at different source temperatures. Values calculated 
at T c = 300K, k L = 0.15 W/m-K. 


power level as a percentage of P max , for example, to 90% of P m ax, 
then for a given Th and ki, the L min to achieve this power output can 
be determined from Eq. (12), a typical result of which is shown in 
Fig. 4(b) for k L = 0.15 W/m-K. 

Figs. 3 and 4 suggest that in order to increase the theoretical 
upper limit of power output at a given T H , L min and k L should be 
taken into consideration during the engine design stage. Taking the 
P3 as an example again, based on the reported dimensions as well 
as Ks and I< L values calculated previously, a theoretical W max of 
9.7 mW can be obtained using Eq. (13). In reality a power output of 
0.35 mW was reported [2], Despite possible errors in the estimation 
of I<s and Ki values, it is still notable that there seems to be plenty of 
room for possible performance improvement in P3 engines. For 
example, the power output of P3 engines have been further 
improved to 2.5 mW with the implementation of active cooling [2], 
Alternatively, improving the thermal efficiency of the thermody¬ 
namic cycle, e.g. by means of a regenerator, could help the engine to 
achieve power generation closer to its thermodynamic maximum. 
In addition, the model can be used to evaluate the relative effects of 
design changes, e.g. adjusting L or ki to reduce leakage, for example 
by increasing the P3's length from 75 pm to 200 pm, the theoretical 
limit of power output predicted by this model would increase from 
9.7 mW to 10.4 mW, a 7% improvement. The relatively small 
improvement of this design change suggests that heat leakage is 
not a major concern in the current deigns of the P3. 

3.2. Power density 

For power applications where mobility, weight, and space are 
major concerns, e.g. as a portable power supply, the power density 
is a key performance objective [lj and thus is discussed in this 
section. The theoretical upper limit of power density is calculated as 
the maximum power divided by the volume, as given in Eq. (15). 
Since weight is generally proportional to volume, the power 
density calculated here is taken to be the equivalent of power per 
unit weight. 

Fig. 5 shows the power density (normalized against the 
maximum power density reachable) which shows that as the 
length scale is reduced, power density increases until it reaches an 
optimum value and then decreases due to the domination of the 
heat leakage. Thus, to aim for high power density with a given 
engine material conductivity ( k L ) and source temperature (Th), it is 
important to design an engine with an optimum length. However, 
in the case of variable Th, the optimum length L op t is a function of 
T h , as shown in Fig. 6 (a). The optimum lengths can be calculated 
from Eq. (14) and a set of typical results is shown in Fig. 6 (b). 


For the case of the P3 engine, the model predicts that the 
theoretical limit of power density is 6.6 mW/mm 3 for the actual 
engine length reported (I = 75 pm). Since the actual density 
reported is only 0.24 mW/mm 3 , as noted in the power discussion, 
further development is needed to bring real power and power 
density closer to these upper limits. Nonetheless, the relative 
effects of a design change can be evaluated by means of the model. 
Based on the reported operating conditions of the P3 (Th = 333K, 
Tc = 300K, and k L = 0.057 W/m-K), the optimal length for 
maximum power density can be calculated to be 4 pm, which 
corresponds to a theoretical limit of maximum power density of 
33.9 mW/mm 3 . We acknowledge that practical fabrication issues 
might limit the realization of such dimensions and power density. 

3.3. Thermal efficiency 

Thermal efficiency for micro heat engines has been the focus of 
studies in previous models. Other factors affect the overall system 
efficiency of micro heat engines, such as mechanical friction and 
fluid-dynamic losses, however only the thermal efficiency is 
considered in this work. 



Fig. 5. Power density vs length scale at different thermal conductivities, power density 
normalized to value without leakage. Values calculated at T H = 400K, T c = 300K. 
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Fig. 6. (a) Power density vs length scale at different source temperatures, power density normalized to value without leakage, (b) Optimum lengths at different source temperature. 
Values calculated at T c = 300K, k L = 0.15 W/m-K. 


Generally speaking, the increased heat leakage due to down- 
scaling of micro heat engines leads to a decrease in thermal effi¬ 
ciency. Although relatively low efficiency micro heat engines can 
still find practical applications [24], it is important to improve 
efficiency especially for power conversion applications [2]. 

The upper limit of thermal efficiency for any heat engine is the 
Carnot efficiency, but it has been shown from heat engine models 
such as the one in Fig. 1 (a) that the thermal efficiency at maximum 
power is somewhat lower than Carnot efficiency [16]. The thermal 
efficiency at max power however has been shown to correspond 
better with real power plant data [14], and hence may be more 
realistic to use for performance evaluation. The effect of length L on 
the thermal efficiency at maximum power ri can be obtained from 
(16), and a set of typical results is shown in Fig. 7. 

For a given Kl, the upper limit of thermal efficiency at maximum 
power can be found at large values of L. Alternatively, since thermal 
efficiency decreases with L, the lower limit of L can be found for 
a targeted efficiency. Note also that choosing a material with lower 
thermal conductivity k L allows a lower value for L, which agrees 
with the result from previous modeling work [11]. For the case of 
the P3 engine, with T H = 333K, Tc = 300K, and k L = 0.057 W/m-K, 
the thermal efficiency at maximum power achievable is 1.47% as 



L (mm) 


Fig. 7. Thermal efficiency vs Length, thermal efficiency normalized against Carnot 
efficiency. Values calculated at T H - 400K, T c = 300K. 


calculated by this model, while the actual thermal efficiency of 0.2% 
is reported [2], Tailoring of k L or L can potentially improve the 
thermal efficiency. For example, if the design length of the P3 were 
changed from 75 pm to 200 pm, the theoretical efficiency can be 
potentially improved from 1.47% to 2.7%. 

3.4. Design trade-off for performance 

Ideally in micro heat engine design, the designed performance 
should include high power, high power density, and high efficiency, 
but in practice micro heat engine performance may involve tradeoffs 
among these objectives. The trade-off between power and thermal 
efficiency is known for general heat engines from endoreversible 
thermodynamics [16]; that is, at maximum Carnot efficiency the 
engine is operating infinitely slow, hence generating zero power. 
Maximum power is generated at a thermal efficiency somewhat 
below the Carnot efficiency. However in real heat engines with 
leakage, it has been shown that the maximum power point and 
maximum thermal efficiency point are in close proximity [15], and 
can practically be treated as equivalent. Hence, maximizing power 
does tend to maximize thermal efficiency, as Eq. (16) shows. 

Another trade-off that needs to be considered for micro heat 
engines is that between power and power density. It is important to 
note that the optimum length for maximum power density is 
different from that for maximum power. The L for maximum power 
is relatively large in order to minimize leakage, however this large L 
(i.e. large volume) leads to reduced power density. Hence, reduced 
power output is obtained when the optimum length for power 
density is designed. One possible way to benefit from high power 
density is to use a large array of micro heat engines with high 
power density to collectively produce high overall power output 
collectively, which would be greater than that of a single engine 
with similar overall volume if losses from interconnection of such 
an array is ignored. Thus, with optimized designs, micro heat 
engine arrays could be collectively used to generate high power 
output if they can be potentially assembled into large areas or 
volumes. 

4. Conclusion 

A new thermodynamic model was developed to evaluate the 
theoretical maximum power output, power density, and thermal 
efficiency of micro heat engines, which includes both heat leakage 
and finite heat input elements. By analyzing this model, it was 
discovered that as the length of heat engines are reduced to 
micron-scale, the theoretical limit of power output generally 
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decreases. In order to improve power output, a lower limit to the 
micro engine length L can be identified for a given set of design 
parameters such as temperature and geometric dimensions; simi¬ 
larly, suitable materials can be selected to reduce detrimental heat 
leakage. Furthermore, the theoretical limit of power density was 
evaluated. For a given engine material (/q.) and source temperature 
(Th), an optimal length can be found that maximizes the limit of 
power density. The theoretical limit of thermal efficiency at 
maximum power for micro heat engines was also studied and 
found to trend similarly to power output. Finally, in aiming for ideal 
overall performance, design trade-offs for power, power density, 
and thermal efficiency were examined, possibly leading to micro 
heat engine arrays that have higher overall output power than 
a single engine with equivalent size. The predictions of the devel¬ 
oped model have been compared to a typical micro heat engine 
throughout the studies of theoretical limits of power, power 
density, and thermal efficiency. By comparing these results, 
potential design improvements were discussed. 
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